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0^ ■ ABSTRACT 



Chandrasekhar separated the Dirac equation for spinning and massive 
particles in Kerr geometry into radial and angular parts. Here we solve the 



> 

04 ' complete wave equation and find out how the Dirac wave scatters off Kerr 
»■ . 

' black holes. The eigenfunctions, eigenvalues and reflection and transmission 

O ! 

, co-efficients are computed. We compare the solutions with several parameters 

. to show how a spinning black hole distinguishes mass and energy of incoming 

I ' waves. Very close to the horizon the solutions become independent of the 

^ particle parameters indicating an universality of the behaviour. 

^ ' Key v^rords: Black holes - spin-1/2 particles - Dirac equations - Waves: 

scattering 



1 INTRODUCTION 

Chandrasekhar (1976) separated Dirac equation in Kerr black hole geometry into radial 
(r) and angular (6) parts. The radial equations governing the radial wave-functions Rj_i 
corresponding to spin ±i are given by (with H = I = G = c): 

A^Po^-i = (A + impr)A^2R^^- A^pjA^i?^i = (A - impr)R_i, (1) 

where, the operators Vn and are given by, 

^ tK (r-M) . ^ iK ir-M) 

V^ = dr + — + 2n^-^^- Vl = dr- — + 2n^-^^, (2) 
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and 

A = + - 2Mr; X = (r^ + a^)(7 + am. (3) 

Here, a is the Kerr parameter, n is an integer, a is the frequency of incident wave, M is the 
mass of the black hole, rrip is the rest mass of the Dirac particle, 7^ is the eigenvalue which 
is the separation constant of complete Dirac equation and m is the azimuthal quantum 
number. 

The equations governing the angular wave-functions Sj_i corresponding to spin ±| are 
given by: 



CiS.i = —{X — ampCos9)S_i; C\S_i = +{X + ampCos9)S,i (4) 

2 2 2 2 ^ ^ 

where, the operators and are given by, 

Cn = de + Q + ncot9; = de - Q + ncotO (5) 

and 

Q — aa sinO + m cosec 9. (6) 

Combining eqs. 4, one obtains a second order angular eigenvalue equations which admits 
exact solutions for spin-half particles when p — ^ — 1 (Chakrabarti, 1984), 



1 - 



2(/ -Fl) aax 



and 



where. 



and 



2(/ + l) + aaa;t 



(7) 



(8) 



p^F{l,iy, x^F{l + 1,1 + 1); y^F{l,l + l) 



F{k,l2) = [(2/2 + l)(2/i + l)]^ < klmOlkm > [< /sl^Ol/i^ > +(-1)^^-^ < klmOlhrn > (9) 

[< kHoiiA > +{-iy^-^pV2 < hi - \i\h\ >]. 



Here, < ....|.. > are the usual Clcbsh-Gordon coefficients and sY/m arc the standard spin- 
weighted spherical harmonics (Chakrabarti, 1984; see also, Goldberg et al 1967, Breure et 
al, 1982) of spin s and usual quantum numbers I and m. When p ^ ^ — 1, one obtains the 
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solutions perturbatively with aa to be the perturbation parameter. The detailed procedure 
to obtain eigenfunctions and eigenvalues is in (Chakrabarti 1984) and is not described here. 

The radial equations (1) are in coupled form. One can decouple them and express the 
equation either in terms of spin up or spin down wave functions Rj_i but the expression 
loses its transparency. It is thus advisable to use the approach of Chandrasekhar (1983) by 

changing the basis and independent variable r to, 

2Mr+ + am/a^ ( r \ 2Mr^ + am/a, ( r \ , , , , 

r, = r + ± ^log 1 ^log 1 r > r+ . 10 

r+ — r_ yr+ j r+ — r_ yr_ j 

where, 

d d 111 22 / 

- — = ^-r] u! = r + a ; a =a +am/a, (11) 
ctr* co^ dr 

to transform the set of coupled equations (eq. 1) into two independent one dimensional wave 
equations given by: 



Here, Vq — ^{-£r and Vl — ^{-£ — icr) were used and wave functions were redefined 



as R_i = P_i and A^R.i = P, i. 

2 2 "^"2 "'"2 



2 SOLUTION PROCEDURE 



We define a new variable, 

9 = tan~^{mpr /X), which gives, (A ± irupr) = exp{±i9)^{X^ + rn^r"^). (13) 
Also define. 



1 

--i tan 
2 



1 ( TT^p'^ 



P-1 — 



1 . _i frripr 
+ -i tan I — — 
2 V A 



(14) 



and further choosing f* = r* + ^tan ^(IMl) so that dr* = (1 + -%4r^ ^2 , r, )dr*, and 
2i± — ip+i ± ■0_ 1 the above equations become. 



dr* 



where. 



A5(A2 + m2r2)3/2 



or* / 



W = 



a;2(A2 + m2r2)+AmpA/2(7 



(15) 



(16) 



^Prom these equations, we readily obtain a pair of independent one-dimensional wave equa- 
tions, 

(17) 



^ + (T^ ) Z± = V±Z±; where V± = W'^± 



dr 



- 2 



dr* 
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By transformation of the variable from r to (and f^,) the horizon is shifted from r = r+ 
to r* = —oo unless cr < (Jg = —am/2Mr^ (eq. 10). In this connection, it is customary to 
define (7c where — (eq. 11). Thus, (Jc = —m/a. If cr < (jg, super-radiation is expected for 
particles with integral spins but not for those with half- integral spin (Chandrasekhar 1983). 
Thus, we concentrate on the region where, a > (Tg- 

The choice of parameters is generally made in such a way that there is a significant 
interaction between the particle and the black hole, i.e., when the Compton wavelength of 
the incoming wave is of the same order as the outer horizon of the Kerr black hole. Similarly, 
the frequency of the incoming particle (or wave) should be of the same order as inverse of 
light crossing time of the radius of the black hole. These yield, 

mp ~ (T ~ [M + ^[{M'^ - a^)]-\ (18) 

Thus, we need to deal with quantum black holes to get 'interesting' results. There are two 
cases of interest: (1) the waves do not 'hit' the potential barrier and (2) the waves do hit 
the potential barrier. First, we replace the potential barrier by a large number of steps as 
in the step-barrier problem in quantum mechanics. Fig. 1 shows one such example of the 
potential barrier V+ (eq. 17) which is drawn for a = 0.5, nip = 0.8 and a = 0.8. In reality 
we use tens of thousands of steps with suitable variable widths so that the steps become 
indistinguishable from the actual function. The solution of (17) at nth step can be written 
as (Davydov 1976), 

= Anexp[iknh^n] + Bnexp[-iknh,n\ (19) 

when energy of the wave is greater than the height of the potential barrier. The standard 
junction condition is given as (Davydov 1976), 

Z+,n^Z+^n+i and = -^U+i- (20) 



The reflection and transmission co-efficients at nth junction are given by: 

r, ^n+l{kn+l ~ kn) + Bn+i{kn+i ^ hn) rr i d rr,l\ 
i<-n—-, 77 , , X , R 77 J^n — i- - i<-n {^i-} 

At each of the n steps these conditions were used to connect solutions at successive steps. 



Here, k is the wave number {k = ^^a"^ — V±) of the wave and A;„ is its value at nth step. We 
use the 'no-reflection' inner boundary condition: ^ at r* — > — oo. 

For the cases where waves hit on the potential barrier, inside the barrier (where cr^ < V^) 
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Fig. 1 : Behaviour of V+ (smooth soUd curve) for a = 0.5, rup = 0.8, (T = 0.8. This is approximated as a collection of steps. 
In realily tens of thousand steps were used with varying step size which mimic the potential with arbitrary accuracy. 

we use the wave function of the form 

Z+^ri = Ancxpl-anr^^n] + 5„exj9[a„f*,„] (22) 
where, = y ~ c"^) as in usual quantum mechanics. 
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Fig. 2a: Reflection (R) and transmission (T) coefliicients of waves with varying mass as functions of ft. nip = 0.78 (solid), 
nip = 0.79 (dotted) and m,p = 0.80 (long-dashed) are used. Other parameters are a = 0.5 and a = 0.8. Inset shows R in 
logarithmic scale which falls off exponentially just outside the horizon. 

3 EXAMPLES OF SOLUTIONS 

Fig. 2a shows three solutions [amphtudes of Re(Z+)] for parameters: a = 0.5, a = 0.8 and 
nip = 0.78, 0.79, and 0.80 respectively in solid, dotted and long-dashed curves. The energy 
0"^ is always higher compared to the height of the potential barrier (Fig. 1) and therefore 
the particles do not 'hit' the barrier, k goes up and therefore the wavelength goes down 
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Fig. 2b: Amplitude of Re(Z+) of waves with varying mass as functions of ft. rup = 0.78 (solid), mp = 0.79 (dotted) and 
rrip = 0.80 (long-dashed) are used. Other parameters are a = 0.5 and a = 0.8. 

monotonically as the wave approaches a black hole. It is to be noted that though ours is 
apparantly a 'crude' method, it has flexibility and is capable of presenting insight into the 
problem, surpassing any other method such as ODE solver packages. This is because one 
can choose (a) variable steps depending on steepness of the potential to ensure uniform 
accuracy, and at the same time (b) virtually infinite number of steps to follow the potential 
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as closely as possible. For instance, in the inset, we show R in logarithmic scale very close 
to the horizon. All the three curves marge, indicating that the solutions are independent 
of the mass of the particle and a closer inspection shows that here, the slope of the curve 
depends only on a. The exponential dependence of i?„ close to the horizon becomes obvious. 
Asymptotically, V± — nip (eq. 17), thus, as rUp goes down, the wavelength goes down. In Fig. 
2b, wc present the instantaneous values of the reflection R and transmission T coefficients 
(i.e., Rn and T„ of eq. 21) for the same three cases. As the particle mass is decreased, k goes 
up and corresponding R goes down consistent with the limit that as k —>■ oo, there would 
be no reflection at all as in a quantum mechanical problem. 

Figs. 3(a-b) compare a few solutions where the incoming particles 'hit' the potential 
barrier. We choose, a — 0.95, a — 0.168 and mass of the particle nip — 0.16, 0.164, 0.168 
respectively in solid, dotted and long-dashed curves. Inside the barrier, the wave decays 
before coming back to a sinusoidal behaviour, before entering into a black hole. In Fig. 
3b, we plotted the potential (shifted by 2.05 along vertical axis for clarity). Here too, the 
reflection coefficient goes down as k goes up consistent with the classical result that as the 
barrier height goes up more and more, reflection is taking place strongly. Note however, that 
the reflection is close to a hundred percent. Tunneling causes only a few percent to be lost 
into the black hole. 

Figs. 4(a-b) show the nature of the complete wave function when both the radial and 
the angular solutions (Chakrabarti 1984) are included. Fig. 4a shows contours of constant 
amplitude of the wave (i?_i/2'S'_i/2) in the meridional plane - X is along radial direction 
in the equatorial plane and Y is along the vertical direction. The parameters are a — 0.5, 
nip — 0.8 and a — 0.8. Some levels are marked. Two successive contours have amplitude 
difference of 0.1. In Fig. 4b a three-dimensional nature of the complete solution is given. 
Both of these figures clearly show how the wavelength varies with distance. Amplitude of 
the spherical wave coming from a large distance also gets weaker along the vertical axis and 
the wave is forced to fall generally along the equatorial plane, possibly due to the dragging 
of the inertial frame. 

4 CONCLUSION 

Scatterings of massive, spin-half particles from a spinning black hole has been studied with 
particular emphasis to the nature of the radial wave functions and the reflection and trans- 
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Fig. 3a: Reflection (R) coefficient of waves with varying mass as functions of r». nip = 0.16 (solid), nip = 0.164 (dotted) and 
rrip = 0.168 (long-dashed) are used. Other parameters are a = 0.95 and <j = 0.168. 

mission coefficients. Well known quantum mechanical step-potential approach is used but we 
applied it successfully to a complex problem of barrier penetration in a spacetime around a 
spinning black hole. Among significant observations, we find that the wave function and R, 
and T behave similarly close to the horizon independent of the initial parameter, such as the 
particle mass rrip. Particles of different mass scatter off to a large distance completely differ- 
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Fig. 3b: Amplitude of Re{Z+) of waves with varying mass as functions of r,. rup = 0.16 (solid), mp = 0.164 (dotted) and 
mp = 0.168 (long-dashed) are used. Nature of potential with nip = 0.168 is drawn shifting vertically by 2.05 unit for clarity. 
Other parameters are a = 0.95 and a = 0.168. 

ently, thus giving an impression that a black hole could be treated as a mass spectrograph! 
When the energy of the particle becomes higher compared to the rest mass, the reflection 
coefficient diminishes as it should it. Similar to a barrier penetration problem, the reflection 
coefficient becomes close to a hundred percent when the wave hits the potential barrier. 
Another significant observation is that the reflection and transmission coefficients are func- 
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Fig. 4a: Contours of constant amplitude arc plotted in the meridional plane around a black hole. Radial direction on equatorial 
plane is along X axis and the vertical direction and along Y. Both radial and theta solutions have been combined. Parameters 
are a = 0.5, rup = 0.8 and cr = 0.8. 

tions of the radial coordinates. This is understood easily because of the very nature of the 
potential barrier which is strongly space dependent which we approximate as a collection 
of steps. Combining with the solution of theta-equation, we find that the wave- amplitude 
vanishes close to the vertical axis, possibly due to the frame-dragging effects. 
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Fig. 4b; Three dimensional view of R-1/2S-1/2 plotted in the meridional plane around a black hole. Both radial and 
theta solutions have been combined. Parameters are a = 0.5, rup = 0.8 and cr = 0.8. 
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